Control method for consensus of agricultural multi-agent system based on sampled data

ABSTRACT

The present invention relates to the field of control engineering, in particular to a distributed control method for consensus of an agricultural multi-agent system based on sampled data, comprising the following steps: for a first-order multi-agent system with fixed directed topology, designing a distributed control protocol based on sampled information in a time delay state; obtaining a dynamic model of the multi-agent system with time delays based on sampled information, and transforming a consensus problem of multiple agents into a stability problem by means of tree transformation; determining constraint conditions on a time delay and a sampling period that the multi-agent system achieves stability, that is, sufficient and necessary conditions that the agents in the multi-agent system achieve average state consensus; and realizing average consensus of the multi-agent system according to the sufficient and necessary conditions.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims foreign priority of Chinese Patent Application No. 202210626977.9, filed on Jun. 6, 2022 in the China National Intellectual Property Administration, the disclosures of all of which are hereby incorporated by reference.

BACKGROUND OF THE INVENTION 1. Technical Field

The present invention relates to the field of control engineering, in particular to a distributed control method for consensus of an agricultural multi-agent system based on sampled data.

2. Description of Related Art

Driven by digital revolution, agriculture has gradually begun to move towards intelligence, and intelligent agriculture is an important development direction of agriculture in the digital age. The intelligent agriculture, which integrates theoretical knowledge of mathematics, automation, computer, information communication and agriculture and applies same to agricultural agents of hardware, has received extensive attention from scholars of related disciplines such as information technology, computer science and agricultural science.

The planting methods of different crops are greatly different, but need a lot of labor, particularly in the aspects of planting, spraying, picking, etc. of crops. Even during farming, the use of chemical pesticides or frequent repetitive actions in a short time may produce a certain impact on the health of workers. Therefore, the level of agricultural mechanization needs to be further improved, and the emergence of multi-agent systems provides a new development trend for reducing labor and improving the level of agricultural mechanization.

A multi-agent system refers to a group system composed of a plurality of independent individuals, and its goal is communication and interaction through mutual information among the individuals. An agent refers to an autonomous individual with the basic characteristics of autonomy, sociality, response and pre-activity, and may be regarded as a corresponding software program or an entity (such as a person, a vehicle, or a robot). Among many researches on multi-agent systems, the consensus problem is one of the most valuable topics. For agricultural multi-agent systems, the consensus problem is more important. In traditional multi-agent system consensus control protocols, a lot of information interaction is required to achieve a state consensus, and the multi-agent systems are required to be equipped with high-performance processors, which improve the design and manufacturing costs of agricultural multi-agent systems.

BRIEF SUMMARY OF THE INVENTION

In order to reduce redundant information interaction and reduce requirements for processor configuration, a distributed control protocol based on sampled data was introduced into the research on the consensus problem of an agricultural multi-agent system. Different from other continuous control protocols, the control protocol based on sampled data reduces information redundancy, reduces the control cost of a system and the requirements of network communication, and improves the robustness and redundancy of the system. Meanwhile, in order to make the control protocol more practical and achieve average state consensus among agents in the case of time delay, the distributed control protocol based on sampled data was considered in the presence of time delay.

The present invention provides the following technical solution: a distributed control method for consensus of an agricultural multi-agent system based on sampled data, including the following steps:

-   -   Step 1, for a first-order multi-agent system with fixed directed         topology, designing a distributed control protocol based on         sampled information in a time delay state;     -   Step 2, obtaining a dynamic model of the multi-agent system with         time delays based on sampled information, and transforming a         consensus problem of multiple agents into a stability problem by         means of tree transformation;     -   Step 3, determining constraint conditions on a time delay and a         sampling period that the multi-agent system achieves stability,         that is, sufficient and necessary conditions that the agents in         the multi-agent system achieve average state consensus; and     -   Step 4, realizing average consensus of the multi-agent system         according to the sufficient and necessary conditions that the         multi-agent system achieves average consensus.

In step 1, for a system including N agents, the states of the agents are represented by x_(i), where i=1, 2, . . . , N; the communication topology directed graph G=(V, E, A) of the networked multi-agent system is a weighted directed graph, N vertices v in the directed graph G represent N agents in the multi-agent system, and the i-th vertex of the directed graph G is represented by v_(i), where i=1, 2, . . . , N; V={v₁, v₂, . . . , v_(N)} represents a set of vertices, and the vertex v_(i) is the i-th vertex in the directed graph G and corresponds to the i-th agent in the multi-agent system; there are totally N vertices, each agent is a vertex of the directed graph G, and the state level of each vertex represents an actual physical value, including position, temperature, or voltage; E⊆V×V is a frontier set, and A=[a_(ij)] is a non-negative weighted adjacency matrix, where j=1, 2, . . . , N; a directed edge from the vertex v_(i) to v_(j) is E_(ij)=(v_(j), v_(i)), an adjacency matrix element a_(ij) with regard to E_(ij) is a non-negative real number, and a set of neighborhood nodes of the vertex v_(i) is N_(i)={v_(i)∈V|(v_(j),v_(i))∈E}; if there is at least one directed edge between two vertices, the directed graph G is a strongly connected graph, and the directed graph G has an indegree matrix

$\Delta = {{diag}\left\{ {{\sum\limits_{j = 1}^{n}a_{1j}},{\sum\limits_{j = 1}^{n}a_{2j}},\ldots,{\sum\limits_{j = 1}^{n}a_{nf}}} \right\}}$

and a Laplacian matrix L=[l_(ij)]∈R^(n×n), where L=Δ−A.

An element l_(ij) in the Laplacian matrix satisfies

$l_{ij} = \left\{ {\begin{matrix} {{\sum\limits_{j \in {N(i)}}a_{ij}},{j = i}} \\ {{- a_{ij}},{j \neq i}} \end{matrix}.} \right.$

Because the agricultural multi-agent system studied is strongly connected, a diagonal matrix W=diag{det(L₁₁), det(L₂₂), . . . , det(L_(nn))} may be derived, a left eigenvector being w_(i)=[det(L₁₁), det(L₂₂), . . . , det(L_(nn))], where ω^(T)L=0_(n) ^(T), L_(ii)∈R^((n-1)×(n-1)) is a matrix after the i-th row and the i-column are removed from the Laplacian matrix, where det(L_(ii)) represents a determinant of the matrix L_(ii); a new topology graph G=(V, Ē, Ā) may be obtained according to the diagonal matrix W, where an element ā_(ij) in Ā=[ā_(ij)] satisfies the following:

${{\overset{\_}{a}}_{ij} = \frac{{\det\left( L_{ii} \right)a_{ij}} + {\det\left( L_{jj} \right)a_{ji}}}{2}},{\forall i},{j \in {\mathcal{I}_{n}.}}$

Where,

_(n)={1, 2, . . . , N}, the following may be further obtained:

$\overset{\_}{A} = {\frac{{WA} + {A^{T}W}}{2}.}$

The relationship between the Laplacian matrix L of the topology graph G and the Laplacian matrix L of the topology graph G is as follows:

$\overset{\_}{L} = {\frac{{WL} + {L^{T}W}}{2}.}$

In addition, all other eigenvalues λ₂, . . . , λ_(n) of the Laplacian matrix L except a zero eigenvalue λ₁ have a positive real part, and the mirror graph G of the directed graph G satisfies the properties of an undirected strongly connected graph.

Considering the average consensus problem of the multi-agent system with strongly connected directed topology, the relationship between agents is represented by an edge relationship between vertices. In the directed graph G of the multi-agent system, the state of each vertex v_(i) is represented by x_(i), a state vector of vertices is represented by x(t), x(t)=[x₁(t), x₂(t), . . . , x_(n)(t)]^(T)∈R^(n), and the dynamic model of the first-order multi-agent system with fixed directed topology is represented as follows:

{dot over (x)} _(i)(t)=u _(i)(t),∀i∈

_(n).

u_(i)(t) is a control input used for solving the consensus problem.

In order to reduce the communication cost of the intelligent agricultural multi-agent system, the objective of the present invention is to solve average consensus problem of the agricultural multi-agent system by using sampled data. In practical applications of the agricultural multi-agent system, the consensus may also be affected by a communication delay. Especially in crop planting, multiple agents need to transmit information to each other in the process of completing a planting task, and excessive communication delay may lead to oscillation or divergence of the multi-agent system, so the problem of time delay needs to be considered. For the communication delay of the agricultural multi-agent system, the sampling period is set as p. Considering the existence of a time delay τ shorter than one sampling period, the proposed distributed delay control protocol based on sampled data is as follows:

${u_{i}(t)} = \left\{ {\begin{matrix} {{\sum\limits_{j \in {N(i)}}{{\overset{\_}{a}}_{ij}\left( {{x_{j}\left( {{kp} - p} \right)} - {x_{i}\left( {{kp} - p} \right)}} \right)}},{t \in \left\lbrack {{kp},{{kp} + \tau}} \right\rbrack}} \\ {{\sum\limits_{j \in {N(i)}}{{\overset{\_}{a}}_{ij}\left( {{x_{j}({kp})} - {x_{i}({kp})}} \right)}},{t \in \left\lbrack {{{kp} + \tau},{{kp} + p}} \right\rbrack}} \end{matrix}.} \right.$ ∀i∈

,k=0,1,2, . . . ,0<τ<p.

In step 2, the specific process of obtaining a dynamic model of the multi-agent system with time delays by means of the distributed control protocol based on sampled data, and transforming a consensus problem of multiple agents into a stability problem by means of tree transformation includes:

-   -   obtaining, according to the proposed distributed delay control         protocol, the dynamic model of the first-order multi-agent         system at the sampling period p and time delay τ:

${\begin{bmatrix} {x\left( {{kp} + p} \right)} \\ {x({kp})} \end{bmatrix} = {\phi\begin{bmatrix} {x({kp})} \\ {x\left( {{kp} - p} \right)} \end{bmatrix}}},{k = 0},1,2,\ldots$ ${{{Where}\phi} = \begin{bmatrix} {I_{n} - {\left( {p - \tau} \right)\overset{\_}{L}}} & {{- \tau}\overset{\_}{L}} \\ I_{n} & 0 \end{bmatrix}},$

Where I is a unit matrix, and I_(n) is an n-order unit matrix.

In order to analyze a convergence problem of the system after the protocol is used, the dynamic model is transformed by means of tree transformation:

y₁(kp) = x₁(kp) y₂(kp) = x₁(kp) − x₂(kp) y₃(kp) = x₁(kp) − x₃(kp) ▫ y_(n)(kp) = x₁(kp) − x_(n)(kp);

An invertible matrix Q is obtained:

$Q = {\begin{bmatrix} 1 & 0 & 0 & \ldots & 0 \\ 1 & {- 1} & 0 & \ldots & 0 \\ 1 & 0 & {- 1} & \ldots & 0 \\  \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & 0 & 0 & \ldots & {- 1} \end{bmatrix} = \begin{bmatrix} {C \in R^{1 \times n}} \\ {E \in R^{{({n - 1})} \times n}} \end{bmatrix}}$

Q⁻¹ is obtained:

$Q^{- 1} = {\begin{bmatrix} 1 & 0 & 0 & \ldots & 0 \\ 1 & {- 1} & 0 & \ldots & 0 \\ 1 & 0 & {- 1} & \ldots & 0 \\  \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & 0 & 0 & \ldots & {- 1} \end{bmatrix} = \text{?}}$ ?indicates text missing or illegible when filed

The following is obtained from Q and Q⁻¹:

${y({kp})}\overset{\Delta}{=}{\begin{bmatrix} {y_{1}({kp})} \\ {y_{2}({kp})} \\  \vdots \\ {y_{n}({kp})} \end{bmatrix} = {{{Qx}({kp})}.}}$

Thus, y(kp+p)=Qx(kp+p), y(kp−p)=Qx(kp−p) and y(kp)=Qx(kp) are obtained, and the following is further obtained:

${\begin{bmatrix} {y\left( {{kp} + p} \right)} \\ {y({kp})} \end{bmatrix} = {\begin{bmatrix} {I_{n} - {\left( {p - \tau} \right)H}} & {{- \tau}H} \\ I_{n} & 0_{n} \end{bmatrix}\begin{bmatrix} {y({kp})} \\ {y\left( {{kp} - p} \right)} \end{bmatrix}}},$ ${{{Where}H} = {{Q\overset{\_}{L}Q^{- 1}} = \begin{bmatrix} 0 & {C\overset{\_}{L}F} \\ 0_{n - 1} & {E\overset{\_}{L}F} \end{bmatrix}}};$

Thus, the system is divided into two subsystems:

$\begin{bmatrix} {y_{1}\left( {{kp} + p} \right)} \\ {y_{1}({kp})} \end{bmatrix} = {{\begin{bmatrix} 1 & 0 \\ 1 & 0 \end{bmatrix}\begin{bmatrix} {y_{1}({kp})} \\ {y_{1}\left( {{kp} - p} \right)} \end{bmatrix}} - {\begin{bmatrix} {\left( {p - \tau} \right)C\overset{\_}{L}F} & {\tau C\overset{\_}{L}F} \\ 0 & 0 \end{bmatrix}\begin{bmatrix} {\hat{y}({kp})} \\ {\hat{y}\left( {{kp} - p} \right)} \end{bmatrix}}}$ and $\begin{bmatrix} {y_{1}\left( {{kp} + p} \right)} \\ {y_{1}({kp})} \end{bmatrix} = {{\begin{bmatrix} {I_{n - 1} - {\left( {p - \tau} \right)E\overset{\_}{L}F}} & {{- \tau}E\overset{\_}{L}F} \\ I_{n - 1} & 0_{n - 1} \end{bmatrix}\begin{bmatrix} {\hat{y}\left( {{kp} + p} \right)} \\ {\hat{y}\left( {{kp} - p} \right)} \end{bmatrix}}.}$ Where, ŷ(kp) = [y₂(kp), y₃(kp), …, y_(n)(kp)]^(T)and ŷ(kp − p) = [y₂(kp − p), y₃(kp − p), …, y_(n)(kp − p)]^(T).

As can be seen, the subsystems after dimension reduction achieve stability, indicating that the whole system achieves consensus.

In step 3, the constraint conditions on the time delay and the sampling period that the multi-agent system achieves stability, that is, the sufficient and necessary conditions that the agents in the multi-agent system achieve average state consensus, are obtained by means of bilinearity and Hurwitz stability criteria, specifically:

The following is obtained by using an invertible matrix T:

${T^{- 1}E\overset{\_}{L}{FT}} = {\Lambda = \begin{bmatrix} \lambda_{2} & * & & \\  & \lambda_{3} & \ddots & \\  & & \ddots & * \\  & & & \lambda_{n} \end{bmatrix}}$

Where λ₂, λ₃, . . . , λ_(n) are non-zero eigenvalues of L and elements, * is 0 or 1, and further,

{tilde over (y)}(kp+p)=T ⁻¹ ŷ(kp+p).

{tilde over (y)}(kp)=T ⁻¹ ŷ(kp)

{tilde over (y)}(kp−p)=T ⁻¹ ŷ(kp−p).

The dimension reduction system is transformed into:

${\begin{bmatrix} {\overset{\sim}{y}\left( {{kp} + p} \right)} \\ {\overset{\sim}{y}({kp})} \end{bmatrix} = {\zeta\begin{bmatrix} {\overset{\sim}{y}\left( {{kp} + p} \right)} \\ {\overset{\sim}{y}\left( {{kp} - p} \right)} \end{bmatrix}}},$ Where ${\zeta\begin{bmatrix} {I_{n - 1} - {\left( {p - \tau} \right)\Lambda}} & {- {\tau\Lambda}} \\ I_{n - 1} & 0_{n - 1} \end{bmatrix}}.$

A characteristic polynomial of is further obtained:

${\det\left( {{sI}_{{2n} - 2} - \zeta} \right)} = {{❘\begin{matrix} {{sI}_{n - 1} - \left\lbrack {I_{n - 1} - {\left( {p - \tau} \right)\Lambda}} \right\rbrack} & {\tau\Lambda} \\ {- I_{n - 1}} & {sI}_{n - 1} \end{matrix}❘} = {{❘\xi ❘}.}}$

The following is obtained from the properties of a block matrix:

${❘\xi ❘} = {{❘{sI}_{n - 1}❘}{❘{{sI}_{n - 1} - \left\lbrack {I_{n - 1} - {\left( {p - \tau} \right)\Lambda}} \right\rbrack + {{I_{n - 1}\left( {sI}_{n - 1} \right)}^{- 1}{\tau\Lambda}{❘{= {{\prod_{i = 2}^{n}\left\lbrack {s^{2} - {\left( {1 - {p\lambda_{1}} + {\tau\lambda}_{1}} \right)s} + {\tau\lambda}_{1}} \right\rbrack} = {\prod_{i = 2}^{n}{{g_{i}(s)}.}}}}}}}}}$

Then, the following can be known by bilinear transformation of

$s = {\frac{z + 1}{z - 1}:}$ f _(i)(z)=pλ _(i) z ²+2(1−τλ_(i))z+(p−2τ)λ_(i)+2.

If f_(i)(z) is Hurwitz-stable and g_(i)(s) is Schur-stable, the stability is determined as follows:

Assuming z=ωι, then:

f _(i)(ω)=−pλ _(i)ω²+2(1−τλ_(i))ωι+(p−2τ)λ_(i)+2,

Its real part is:

f _(ω)(ω)=−pλ _(i)ω²+(p−2τ)λ_(i)+2 and

Its imaginary part is:

f _(i)(ω)=2(1−τλ_(i))ω;

Upon verification, f_(i)(z) is Hurwitz-stable when satisfying the following conditions:

${\tau < \frac{1}{\lambda_{\max}}},{0 < p < {\frac{2}{\lambda_{\max}} + {2\tau}}},$

That is, under the Hurwitz-stable condition, the multi-agent system achieves consensus.

Step 4 includes the following steps: because L has an eigenvalue of 0, ϕ has corresponding determinable eigenvalues of 0 and 1; the condition of Hurwitz stability of the system is determined after the tree transformation; when the constraint conditions are satisfied, the modulus values of remaining unknown eigenvalues of ϕ are less than 1, that is, the eigenvalues of ϕ except 0 and 1 are within a unit circle under the conditions; because G has undirected strongly connected, so a right eigenvector ω_(r) and a left eigenvector ω_(l) of ϕ when the eigenvalue is 1 satisfy the following:

${w_{r} = {w_{l} = {\frac{1}{\sqrt{n}}1_{n}}}},$

-   -   and ω_(l) ^(T) ω_(r)=1 is satisfied;

All other eigenvalues of ϕ are within the unit circle and there is:

${\begin{bmatrix} {x\left( {{kp} + p} \right)} \\ {x({kp})} \end{bmatrix} = {{\phi\begin{bmatrix} {x({kp})} \\ {x\left( {{kp} - p} \right)} \end{bmatrix}} = {\phi^{x}\begin{bmatrix} {x(0)} \\ {x(0)} \end{bmatrix}}}},$ So: ${\lim\limits_{k\rightarrow\infty}\begin{bmatrix} {x\left( {{kp} + p} \right)} \\ {x({kp})} \end{bmatrix}} = {{\text{?}\begin{bmatrix} {x(0)} \\ {x(0)} \end{bmatrix}} = \text{?}}$ ?indicates text missing or illegible when filed

Therefore, through the control protocol based on sampled information proposed by the present invention, the agricultural multi-agent system with a first-order dynamic model can also achieve average consensus despite a time delay.

It can be seen from the above description that the average consensus distributed control protocol for the agricultural multi-agent system with time delays based on sampled data in the present invention considers the influence of the time continuous control protocol and time delay on the performance of the multi-agent system. The control protocol is designed by a sampling means and a sampling period, thus obtaining an upper limit of sampling period to achieve average consensus of the system. Compared with a continuous control protocol, the robustness and information utilization of the multi-agent system are improved, and the requirements for system hardware are reduced. The control protocol proposed by the present invention overcomes the influence of time delay on the average consensus. With the help of analysis methods of graph theory and matrix theory, sufficient and necessary conditions that the agricultural multi-agent system achieves average consensus are given, and a sampling period and upper and lower limits of upper and lower limits of a time delay are obtained. Because the protocol in the present invention can enable the multi-agent system to achieve the average consensus, the state value for finally achieving the average consensus can be further adjusted by setting the initial state value of each agent in the agricultural system.

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

FIG. 1 is a flowchart of the present invention.

FIG. 2 is a communication topology diagram of a multi-agent system in a specific embodiment.

FIG. 3 shows a convergence state of each agent when the time delay τ=0.2 and the sampling period p=0.93 in the specific embodiment.

FIG. 4 shows a convergence state of each agent when the time delay τ=0.2616 and the sampling period p=1.04 in the specific embodiment.

FIG. 5 shows a convergence state of each agent when the time delay τ=0.2 and the sampling period p=0.92 in the specific embodiment.

DETAILED DESCRIPTION OF THE INVENTION

Technical solutions in specific embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the specific embodiments of the present invention. Apparently, the specific embodiments described are only some of the specific embodiments of the present invention, not all of them. Based on the specific embodiments in the present invention, all other specific embodiments obtained by those of ordinary skill in the art without any creative effort fall within the protection scope of the present invention.

As can be seen from the accompanying drawings, a distributed control method for consensus of an agricultural multi-agent system based on sampled data in the present invention includes the following four steps.

Step 1, for a first-order multi-agent system with fixed directed topology, a distributed control protocol based on sampled information in a time delay state is designed.

In this step, for a system including N agents, the states of the agents are represented by x_(i), where i=1, 2, . . . , N; the communication topology directed graph G=(V, E, A) of the networked multi-agent system is a weighted directed graph, N vertices v in the directed graph G represent N agents in the multi-agent system, and the i-th vertex of the directed graph G is represented by v_(i), where i=1, 2, . . . , N; V={v₁, v₂, . . . , v_(N)} represents a set of vertices, and the vertex v_(i) is the i-th vertex in the directed graph G and corresponds to the i-th agent in the multi-agent system; there are totally N vertices, each agent is a vertex of the directed graph G, and the state level of each vertex represents an actual physical value, including position, temperature, or voltage; E⊆V×V is a frontier set, and A=[α_(ij)] is a non-negative weighted adjacency matrix, where j=1, 2, . . . , N; a directed edge from the vertex v_(i) to v_(j) is E_(ij)=(v_(j), v_(i)), an adjacency matrix element α_(ij) with regard to E_(ij) is a non-negative real number, and a set of neighborhood nodes of the vertex v_(i) is N_(i)={v_(i)∈V|(v_(j),v_(i))∈E}; if there is at least one directed edge between two vertices, the directed graph G is a strongly connected graph, and the directed graph G has an indegree matrix

$\Delta = {{diag}\left\{ {{\sum\limits_{j = 1}^{n}a_{1j}},{\sum\limits_{j = 1}^{n}a_{2j}},\ldots,{\sum\limits_{j = 1}^{n}a_{nj}}} \right\}}$

and a Laplacian matrix L=[l_(ij)]∈R^(n×n), where L=Δ−A.

An element l_(ij) in the Laplacian matrix satisfies

$I_{ij} = \left\{ {\begin{matrix} {{\sum\limits_{j \in {N(i)}}a_{ij}},{j = i}} \\ {{- a_{ij}},{j \neq i}} \end{matrix}.} \right.$

Because the agricultural multi-agent system studied is strongly connected, a diagonal matrix W=diag{det(L₁₁), det(L₂₂), . . . , det(L_(nn))} may be derived, a left eigenvector being w_(i)=[det(L₁₁), det(L₂₂) . . . , det(L_(nn))], where ω^(T) L=0_(n) ^(T); L_(ii)∈R^((n-1)×(n-1)) is a matrix after the i-th row and the i-column are removed from the Laplacian matrix, where det(L_(ii)) represents a determinant of the matrix L_(ii); a new topology graph G=(V, Ē, Ā) may be obtained according to the diagonal matrix W, where an element in ā_(ij) in Ā=[ā_(ij)] satisfies the following:

${{\overset{\_}{a}}_{ij} = \frac{{{\det\left( L_{ii} \right)}a_{ij}} + {{\det\left( L_{jj} \right)}a_{ji}}}{2}},{\forall i},{j \in \mathcal{I}_{n}}$

Where

_(n)={1, 2, . . . , N}, the following may be further obtained:

$\overset{\_}{A} = {\frac{{WA} + {A^{T}W}}{2}.}$

The relationship between the Laplacian matrix L of the topology graph G and the Laplacian matrix L of the topology graph G is as follows:

$\overset{\_}{L} = {\frac{{WA} + {L^{T}W}}{2}.}$

In addition, all other eigenvalues λ₂, . . . , λ_(n) of the Laplacian matrix L except a zero eigenvalue λ₁ have a positive real part, and the mirror graph G of the directed graph G satisfies the properties of an undirected strongly connected graph.

Considering the average consensus problem of the multi-agent system with strongly connected directed topology, the relationship between agents is represented by an edge relationship between vertices. In the directed graph G of the multi-agent system, the state of each vertex v_(i) is represented by x_(i), a state vector of vertices is represented by x(t), x(t)=[x₁(t), x₂(t), . . . , x_(n)(t)]^(T)∈R^(n), and a dynamic model of the first-order multi-agent system with fixed directed topology is represented as follows:

{dot over (x)} _(i)(t)=u _(i)(t),∀i∈

_(n).

u_(i)(t) is a control input used for solving the consensus problem.

In order to reduce the communication cost of the intelligent agricultural multi-agent system, the objective of the present invention is to solve average consensus problem of the agricultural multi-agent system by using sampled data. In practical applications of the agricultural multi-agent system, the consensus may also be affected by a communication delay. Especially in crop planting, multiple agents need to transmit information to each other in the process of completing a planting task, and excessive communication delay may lead to oscillation or divergence of the multi-agent system, so the problem of time delay needs to be considered. For the communication delay of the agricultural multi-agent system, a sampling period is set as p. Considering the existence of a time delay τ shorter than one sampling period, the proposed distributed delay control protocol based on sampled data is as follows:

${u_{i}(t)} = \left\{ \begin{matrix} {{\sum\limits_{j \in {N(i)}}{{\overset{\_}{a}}_{ij}\left( {{x_{j}\left( {{kp} - p} \right)} - {x_{i}\left( {{kp} - p} \right)}} \right)}},} & {t \in \left\lbrack {{kp},{{kp} + \tau}} \right\rbrack} \\ {{\sum\limits_{j \in {N(i)}}{{\overset{\_}{a}}_{ij}\left( {{x_{j}({kp})} - {x_{i}({kp})}} \right)}},} & {t \in \left\lbrack {{{kp} + \tau},{{kp} + p}} \right\rbrack} \end{matrix} \right.$ ∀i ∈ ℐ_(n), k = 0, 1, 2, …, 0 < τ < p.

Step 2, a dynamic model of the multi-agent system with time delays based on sampled information is obtained, and a consensus problem of multiple agents is transformed into a stability problem by means of tree transformation.

In step 2, the specific process of obtaining a dynamic model of the multi-agent system with time delays by means of the distributed control protocol based on sampled data, and transforming a consensus problem of multiple agents into a stability problem by means of tree transformation includes:

According to the proposed distributed delay control protocol, the dynamic model of the first-order multi-agent system at the sampling period p and time delay τ is obtained:

${\begin{bmatrix} {x\left( {{kp} + p} \right)} \\ {x({kp})} \end{bmatrix} = {\phi\begin{bmatrix} {x({kp})} \\ {x\left( {{kp} - p} \right)} \end{bmatrix}}},{k = 0},1,2,\ldots$ ${{Where}\phi} = {\begin{bmatrix} {I_{n} - {\left( {p - \tau} \right)\overset{\_}{L}}} & {{- \tau}\overset{\_}{L}} \\ I_{n} & 0 \end{bmatrix}.}$

Where I is a unit matrix, and I_(n) is an n-order unit matrix.

In order to analyze a convergence problem of the system after the protocol is used, the dynamic model is transformed by means of tree transformation:

y ₁(kp)=x ₁(kp)

y ₂(kp)=x ₁(kp)−x ₂(kp)

y ₃(kp)=x ₁(kp)−x ₃(kp)

□

y _(n)(kp)=x ₁(kp)−x _(n)(kp)

An invertible matrix Q is obtained:

$Q = {\begin{bmatrix} 1 & 0 & 0 & \ldots & 0 \\ 1 & {- 1} & 0 & \ldots & 0 \\ 1 & 0 & {- 1} & \ldots & 0 \\  \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & 0 & 0 & \ldots & {- 1} \end{bmatrix} = \begin{bmatrix} {C \in R^{1 \times n}} \\ {E \in R^{{({n - 1})} \times n}} \end{bmatrix}}$

Q⁻¹ is obtained:

$Q^{- 1} = {\begin{bmatrix} 1 & 0 & 0 & \ldots & 0 \\ 1 & {- 1} & 0 & \ldots & 0 \\ 1 & 0 & {- 1} & \ldots & 0 \\  \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & 0 & 0 & \ldots & {- 1} \end{bmatrix} = \left\lbrack {w_{r} \in {R^{n \times 1}F} \in R^{n \times {({n - 1})}}} \right\rbrack}$

The following is obtained from Q and Q⁻¹:

${y({kp})}\overset{\Delta}{=}{\begin{bmatrix} {y_{1}(p)} \\ {y_{2}({kp})} \\  \vdots \\ {y_{ii}({kp})} \end{bmatrix} = {{{Qx}({kp})}.}}$

Thus, y(kp+p)=Qx(kp+p), y(kp−p)=Qx(kp−p) and y(kp)=Qx(kp) are obtained, and the following is further obtained:

${\begin{bmatrix} {y\left( {{kp} + p} \right)} \\ {y({kp})} \end{bmatrix} = {\begin{bmatrix} {I_{n} - {\left( {p - \tau} \right)H}} & {{- \tau}H} \\ I_{n} & 0_{n} \end{bmatrix}\begin{bmatrix} {y({kp})} \\ {y\left( {{kp} - p} \right)} \end{bmatrix}}},$ ${{{Where}H} = {Q\overset{\_}{L}{Q^{- 1}\begin{bmatrix} 0 & {C\overset{\_}{L}F} \\ 0_{n - 1} & {E\overset{\_}{L}F} \end{bmatrix}}}};$

Thus, the system is divided into two subsystems:

$\begin{bmatrix} {y\left( {{kp} + p} \right)} \\ {y({kp})} \end{bmatrix} = {{\begin{bmatrix} 1 & 0 \\ 1 & 0 \end{bmatrix}\begin{bmatrix} {y_{1}({kp})} \\ {y_{1}\left( {{kp} - p} \right)} \end{bmatrix}} - {\begin{bmatrix} {\left( {p - \tau} \right)C\overset{\_}{L}F} & {{\tau C}\overset{\_}{L}F} \\ 0 & 0 \end{bmatrix}\begin{bmatrix} {\hat{y}({kp})} \\ {\hat{y}\left( {{kp} - p} \right)} \end{bmatrix}}}$ ${{and}\begin{bmatrix} {y_{1}\left( {{kp} + p} \right)} \\ {y_{1}({kp})} \end{bmatrix}} = {{\begin{bmatrix} {I_{n - 1} - {\left( {p - \tau} \right)E\overset{\_}{L}F}} & {{- \tau}E\overset{\_}{L}F} \\ I_{n - 1} & 0_{n - 1} \end{bmatrix}\begin{bmatrix} {\hat{y}({kp})} \\ {\hat{y}\left( {{kp} - p} \right)} \end{bmatrix}}.}$

Where,

ŷ(kp)=[y ₂(kp),y ₃(kp), . . . ,y _(n)(kp)]^(T) and

ŷ(kp−p)=[y ₂(kp−p),y ₃(kp−p), . . . ,y _(n)(kp−p)]^(T).

As can be seen, the subsystems after dimension reduction achieve stability, indicating that the whole system achieves consensus.

Step 3, constraint conditions on a time delay and a sampling period that the multi-agent system achieves stability, that is, sufficient and necessary conditions that the agents in the multi-agent system achieve average state consensus, are determined.

In step 3, the constraint conditions on the time delay and the sampling period that the multi-agent system achieves stability, that is, the sufficient and necessary conditions that the agents in the multi-agent system achieve average state consensus, are obtained by means of bilinearity and Hurwitz stability criteria, specifically:

The following is obtained by using an invertible matrix T:

${T^{- 1}E\overset{\_}{L}{FT}} = {\Lambda = \begin{bmatrix} \lambda_{2} & * & & \\  & \lambda_{3} & \ddots & \\  & & \ddots & * \\  & & & \lambda_{n} \end{bmatrix}}$

Where λ₂, λ₃, . . . , λ_(n) are non-zero eigenvalues of L and elements, * is 0 or 1, and further,

{tilde over (y)}(kp+p)=T ⁻¹ ŷ(kp+p)

{tilde over (y)}(kp)=T ⁻¹ ŷ(kp)

{tilde over (y)}(kp−p)=T ⁻¹ ŷ(kp−p).

The dimension reduction system is transformed into:

${\begin{bmatrix} {\overset{\sim}{y}\left( {{kp} + p} \right)} \\ {\overset{\sim}{y}({kp})} \end{bmatrix} = {\zeta\begin{bmatrix} {\overset{\sim}{y}\left( {{kp} + p} \right)} \\ {\overset{\sim}{y}\left( {{kp} - p} \right)} \end{bmatrix}}},$ ${{Where}\zeta} = {\begin{bmatrix} {I_{n - 1} - {\left( {p - \tau} \right)\Lambda}} & {{- \tau}\Lambda} \\ I_{n - 1} & 0_{n - 1} \end{bmatrix}.}$

A characteristic polynomial of is further obtained:

${\det\left( {{sI}_{{2n} - 2} - \zeta} \right)} = {{❘\begin{matrix} {{sI}_{n - 1} - \left\lbrack {I_{n - 1} - {\left( {p - \tau} \right)\Lambda}} \right\rbrack} & {\tau\Lambda} \\ {- I_{n - 1}} & {sI}_{n - 1} \end{matrix}❘} = {{❘\xi ❘}.}}$

The following is obtained from the properties of a block matrix:

$\begin{matrix} {{❘\xi ❘} = {{❘{sI}_{n - 1}❘}{❘{{sI}_{n - 1} - \left\lbrack {I_{n - 1} - {\left( {p - \tau} \right)\Lambda}} \right\rbrack + {{I_{n - 1}\left( {sI}_{n - 1} \right)}^{- 1}{\tau\Lambda}}}❘}}} \\ \left. {= {\prod_{i = 1}^{n}\left\lbrack {s^{2} - {\left( {1 - {p\lambda_{l}} + {\tau\lambda_{l}}} \right)s} + {\lambda\Lambda_{l}}} \right.}} \right) \\ {= {\prod_{i = 2}^{n}{{g_{i}(s)}.}}} \end{matrix}$

Then, the following can be known by bilinear transformation of

$s = {\frac{z + 1}{z - 1}:}$ f _(i)(z)=pλ _(i) z ²+2(1−τλ_(i))z+(p−2τ)λ_(i)+2.

If f_(i)(z) is Hurwitz-stable and g_(i)(s) is Schur-stable, the stability is determined as follows:

Assuming z=ωι, then:

f _(i)(ω)=−pλ _(i)ω²+2(1−τλ_(i))ωι+(p−2τ)λ_(i)+2,

Its real part is:

f _(ω)(ω)=−pλ _(i)ω²+(p−2τ)λ_(i)+2 and

Its imaginary part is:

f _(i)(ω)=2(1−τλ_(i))ω;

Upon verification, f_(i)(z) is Hurwitz-stable when satisfying the following conditions:

${\tau < \frac{1}{\lambda_{\max}}},{0 < p < {\frac{2}{\lambda_{\max}} + {2\tau}}},$

That is, under the Hurwitz-stable condition, the multi-agent system achieves consensus.

Step 4, average consensus of the multi-agent system is realized according to the sufficient and necessary conditions that the multi-agent system achieves average consensus.

Step 4 includes the following steps: because L has an eigenvalue of 0, ϕ has corresponding determinable eigenvalues of 0 and 1; the condition of Hurwitz stability of the system is determined after the tree transformation; when the constraint conditions are satisfied, the modulus values of remaining unknown eigenvalues of ϕ are less than 1, that is, the eigenvalues of ϕ except 0 and 1 are within a unit circle under the conditions; because G is undirected strongly connected, so a right eigenvector ω_(r) and a left eigenvector ω_(l) of ϕ when the eigenvalue is 1 satisfy the following:

${w_{r} = {w_{l} = {\frac{1}{\sqrt{n}}1_{n}}}},$

and ω_(l) ^(T) ω_(r)=1 is satisfied;

All other eigenvalues of ϕ are within the unit circle and there is:

${\begin{bmatrix} {x\left( {{kp} + p} \right)} \\ {x({kp})} \end{bmatrix} = {{\phi\begin{bmatrix} {x({kp})} \\ {x\left( {{kp} - p} \right)} \end{bmatrix}} = {\phi^{k}\begin{bmatrix} {x(0)} \\ {x(0)} \end{bmatrix}}}},$ ${{So}:{\lim\limits_{k\rightarrow\infty}\begin{bmatrix} {x\left( {{kp} + p} \right)} \\ {x({kp})} \end{bmatrix}}} = {{\begin{bmatrix} {w_{r}w_{l}^{T}} & 0 \\ {w_{r}w_{l}^{T}} & 0 \end{bmatrix}\begin{bmatrix} {x(0)} \\ {x(0)} \end{bmatrix}} = {\begin{bmatrix} {\frac{1}{n}{\sum}_{i = 1}^{n}{x_{i}(0)}} \\ {\frac{1}{n}{\sum}_{i = 1}^{n}{x_{i}(0)}} \end{bmatrix}.}}$

Therefore, through the control protocol based on sampled information proposed by the present invention, the agricultural multi-agent system with a first-order dynamic model can also achieve average consensus despite a time delay.

When this method is implemented, the communication topology of the agricultural first-order multi-agent system including six agents is represented by G₁, and the states of the six agents are represented by x_(i), x₂, x₃, x₄, x₅, and x₆. It can be seen that the graph G₁ is directed unbalanced, and the weight of edges is 1, where V={v₁, v₂, v₃, v₄, v₅, v₆}, and x=[x₁, x₂, x₃, x₄, x₅, x₆]^(T). The communication topology is shown in FIG. 2 .

Parameters of the system are as follows:

$A = \begin{bmatrix} 0 & 0 & 0 & 1 & 0 & 1 \\ 1 & 0 & 0 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 & 1 & 0 \end{bmatrix}$

A Laplacian matrix of the system is as follows:

$L = \begin{bmatrix} 2 & 0 & 0 & {- 1} & 0 & {- 1} \\ {- 1} & 1 & 0 & 0 & 0 & 0 \\ 0 & {- 1} & 1 & 0 & 0 & 0 \\ 0 & 0 & {- 1} & 1 & 0 & 0 \\ 0 & 0 & 0 & {- 1} & 1 & 0 \\ 0 & 0 & 0 & 0 & {- 1} & 1 \end{bmatrix}$

An initial state value of the system is as follows:

x(0)=[3,4,7,6,5,5]^(T).

The following is obtained by using the principle of a mirror graph:

$L = {\begin{bmatrix} 2 & {- 1} & 0 & {- 0.5} & 0 & {- 0.5} \\ {- 1} & 2 & {- 1} & 0 & 0 & 0 \\ 0 & {- 1} & 2 & {- 1} & 0 & 0 \\ 0 & 0 & {- 1} & 2 & {- 0.5} & {- 0.5} \\ 0 & 0 & 0 & {- 0.5} & 1 & {- 0.5} \\ {- 0.5} & 0 & 0 & 0 & {- 0.5} & 1 \end{bmatrix}.}$

The solution is performed by means of the control protocol of the present invention, and finally sufficient and necessary conditions that the system can achieve average consensus are obtained:

${\tau < \frac{1}{\lambda_{\max}}},{0 < p < {\frac{2}{\lambda_{\max}} + {2\tau}}}$

Where λ_(i) is an eigenvalue of the Laplacian matrix L. It can be seen that the value of time delay τ depends on the eigenvalue of the Laplacian matrix L of the mirror graph, and the value of sampling period p is not only related to the eigenvalue of the Laplacian matrix L of the mirror graph, but also related to the value of the time delay. Because the maximum eigenvalue λ_(max) of the Laplacian matrix L is 3.8229, the upper limit of time delay is calculated to be 0.2616 according to the sufficient and necessary conditions to achieve average consensus, and the upper limit of sampling period p is further calculated according to the actual value of time delay τ.

FIG. 3 is an example where the time delay constraint condition is satisfied, but the period constraint condition is not satisfied. If the value of time delay τ is 0.2, the upper limit of sampling period p is 0.9232 by further calculation according to

$p < {\frac{2}{\lambda_{\max}} + {2{\tau.}}}$

The value of p is 0.93, which exceeds the upper limit. The convergence state value of each agent is shown in FIG. 3 . It can be seen that if the time delay satisfies the constraint condition, but the sampling period does not satisfy the period constraint condition, the agricultural first-order multi-agent system cannot achieve average consensus.

FIG. 4 is an example where the time delay constraint condition is not satisfied, but the period constraint condition corresponding to the time delay is satisfied. It can be known according to the sufficient and necessary conditions of average consensus obtained by the present invention that the upper limit of time delay τ is 0.2616. The current value of time delay is 0.262, which exceeds the upper limit. By further calculation, the upper limit of sampling period p is 1.0472, and the value of p is 1.04. The convergence state value of each agent is shown in FIG. 4 . It can be seen that if the time delay does not satisfy the constraint condition, but the sampling period satisfies the period constraint condition, the agricultural first-order multi-agent system cannot achieve average consensus.

FIG. 5 is an example where both the time delay constraint condition and the period constraint condition are satisfied. The value of time delay is 0.2, which is less than the upper limit; and the value of sampling period p is 0.92, which is less than the upper limit 0.9232 of sampling period when the time delay is 0.2. The convergence state value of each agent is shown in FIG. 5 . It can be seen that if the time delay satisfies the constraint condition and the sampling period also satisfies the period constraint condition, that is, when the values of τ and p the sufficient and necessary conditions of the control protocol of the present invention, the agricultural first-order multi-agent system can achieve average consensus.

As can be seen from FIG. 3 to FIG. 5 , for the agricultural multi-agent system with time delays, the invented distributed control protocol for the average consensus of the agricultural multi-agent system with time delays based on sampled data can effectively achieve average consensus of the system.

Although the specific embodiments of the present invention have been shown and described, it could be appreciated by those of ordinary kill in the art that many changes, modification, substitutions and variations may be made to these embodiments without departing from the principle and spirit of the present invention, and the scope of the present invention is defined by the appended claims and equivalents thereof. 

What is claimed is:
 1. A distributed control method for consensus of an agricultural multi-agent system based on sampled data, comprising the following steps: step 1, designing a distributed control protocol based on sampled information in a time delay state, for a first-order multi-agent system under fixed directed topology; step 2, obtaining a dynamic model of the multi-agent system with time delays based on sampled information, and transforming a consensus problem of multiple agents into a stability problem by means of tree transformation; step 3, determining constraint conditions on a time delay and a sampling period that the multi-agent system achieves stability, that is, sufficient and necessary conditions that the agents in the multi-agent system achieve average state consensus; and step 4, realizing average consensus of the multi-agent system according to the sufficient and necessary conditions that the multi-agent system achieves average consensus; wherein in step 1, for a system including N agents, states of the agents are represented by x_(i), where i=1, 2, . . . , N; a communication topology directed graph G=(V, E, A) of a networked multi-agent system is a weighted directed graph, N vertices v in the directed graph G represent N agents in the multi-agent system, and a i-th vertex of the directed graph G is represented by v_(i), where i=1, 2, . . . , N; V={v₁, v₂, . . . , v_(N)} represents a set of vertices, and the vertex v_(i) is the i-th vertex in the directed graph G and corresponds to a i-th agent in the multi-agent system; there are totally N vertices, each agent is a vertex of the directed graph G, and a state level of each vertex represents an actual physical value, including position, temperature, or voltage; E⊆V×V is a frontier set, and A=[a_(ij)] is a non-negative weighted adjacency matrix, where j=1, 2, . . . , N; a directed edge from the vertex v_(i) to v_(j) is E_(ij)=(v_(j), v_(i)), an adjacency matrix element a_(ij) with regard to E_(ij) is a non-negative real number, and a set of neighborhood nodes of the vertex v_(i) is N_(i)={v_(i)∈V|(v_(j),v_(i))∈E}; when there is at least one directed edge between two vertices, the directed graph G is set as a strongly connected graph, and the directed graph G has an indegree matrix $\Delta = {{diag}\left\{ {{\sum\limits_{j = 1}^{n}a_{1j}},{\sum\limits_{j = 1}^{n}a_{2j}},\ldots,{\sum\limits_{j = 1}^{n}a_{nj}}} \right\}}$ and a Laplacian matrix L=[l_(ij)]∈R^(n×n), where L=Δ−A; an element l_(ij) in the Laplacian matrix satisfies $l_{ij} = \left\{ {\begin{matrix} {{\sum\limits_{j \in {N(i)}}a_{ij}},{j = i}} \\ {{- a_{ij}},{j \neq i}} \end{matrix};} \right.$ for a strongly connected agricultural multi-agent system, a diagonal matrix is W=diag{det(L₁₁), det(L₂₂), . . . , det(L_(nn))}, and a left eigenvector is w_(i)=[det(L₁₁) det(L₂₂), . . . , det(L_(nn))], where ω^(T) L=0_(n) ^(T), L_(ii)∈R^((n-1)(n-1)) is a matrix after the i-th row and the i-column are removed from the Laplacian matrix, where det(L_(ii)) represents a determinant of the matrix L_(ii); a new topology graph G=(V, Ē, Ā) may be obtained according to the diagonal matrix W, where an element ā_(ij) in Ā=[ā_(ij)] satisfies the following: ${{\overset{\_}{a}}_{ij} = \frac{{{\det\left( L_{ii} \right)}a_{ij}} + {{\det\left( L_{jj} \right)}a_{ji}}}{2}},{\forall i},{{j \in \mathcal{I}_{n}};}$ where

_(n)={1, 2, . . . , N}, the following may be obtained after putting [ā_(ij)] into: ${\overset{\_}{A} = \frac{{WA} + {A^{T}W}}{2}};$ a relationship between the Laplacian matrix L of the topology graph G and the Laplacian matrix L of the topology graph G is as follows: ${\overset{\_}{L} = \frac{{WA} + {L^{T}W}}{2}};$ in the directed graph G of the multi-agent system, the state of each vertex v_(i) is represented by x_(i), a state vector of vertices is represented by x(t)=[x₁(t), x₂(t), . . . , x_(n)(t)]^(T)∈R^(n), and the dynamic model of the first-order multi-agent system with fixed directed topology is represented as follows: {dot over (x)} _(i)(t)=u _(i)(t),∀i∈

_(n), where u_(i)(t) is a control input used for solving a consensus problem; for a communication delay of the agricultural multi-agent system, a sampling period is set as p; considering an existence of a time delay τ shorter than one sampling period, a proposed distributed delay control protocol based on sampled data is as follows: ${u_{i}(t)} = \left\{ {\begin{matrix} {{\sum\limits_{j \in {N(i)}}{{\overset{\_}{a}}_{ij}\left( {{x_{j}\left( {{kp} - p} \right)} - {x_{i}\left( {{kp} - p} \right)}} \right)}},{t \in \left\lbrack {{kp},{{kp} + \tau}} \right\rbrack}} \\ {{\sum\limits_{j \in {N(i)}}{{\overset{\_}{a}}_{ij}\left( {{x_{j}({kp})} - {x_{i}({kp})}} \right)}},{t \in \left\lbrack {{{kp} + \tau},{{kp} + p}} \right\rbrack}} \end{matrix},{where}} \right.$ ∀i ∈ ℐ_(n), k = 0, 1, 2, …, 0 < τ < p.
 2. The distributed control method for consensus of an agricultural multi-agent system based on sampled data according to claim 1, wherein in step 2, a specific process of obtaining a dynamic model of the multi-agent system with time delays by means of the distributed control protocol based on sampled data, and transforming a consensus problem of multiple agents into a stability problem by means of tree transformation comprises: obtaining, according to a proposed distributed delay control protocol, the dynamic model of the first-order multi-agent system at the sampling period p and time delay τ: ${\begin{bmatrix} {x\left( {{kp} + p} \right)} \\ {x({kp})} \end{bmatrix} = {\phi\begin{bmatrix} {x({kp})} \\ {x\left( {{kp} - p} \right)} \end{bmatrix}}},$ k = 0, 1, 2, … where ${\phi = \begin{bmatrix} {I_{n} - {\left( {p - \tau} \right)\overset{\_}{L}}} & {{- \tau}\overset{\_}{L}} \\ I_{n} & 0 \end{bmatrix}},$ where I is a unit matrix, and I_(n) is an n-order unit matrix.
 3. The distributed control method for consensus of an agricultural multi-agent system based on sampled data according to claim 2, wherein in order to analyze a convergence problem of the system after the protocol is used, the dynamic model is transformed by means of tree transformation: $\begin{matrix} {{y_{1}({kp})} = {x_{1}({kp})}} \\ {{y_{2}({kp})} = {{x_{1}({kp})} - {x_{2}({kp})}}} \\ {{y_{3}({kp})} = {{x_{1}({kp})} - {x_{3}({kp})}}} \\  \vdots \\ {{{y_{n}({kp})} = {{x_{1}({kp})} - {x_{n}({xp})}}};} \end{matrix}$ an invertible matrix Q is obtained: ${Q = {\begin{bmatrix} 1 & 0 & 0 & \cdots & 0 \\ 1 & {- 1} & 0 & \cdots & 0 \\ 1 & 0 & {- 1} & \cdots & 0 \\  \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & 0 & 0 & \cdots & {- 1} \end{bmatrix} = \begin{bmatrix} {C \in R^{1 \times n}} \\ {E \in R^{{({n - 1})} \times n}} \end{bmatrix}}};$ Q⁻¹ is obtained: ${Q^{- 1} = {\begin{bmatrix} 1 & 0 & 0 & \cdots & 0 \\ 1 & {- 1} & 0 & \cdots & 0 \\ 1 & 0 & {- 1} & \cdots & 0 \\  \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & 0 & 0 & \cdots & {- 1} \end{bmatrix} = \begin{bmatrix} {w,{\in R^{n \times 1}}} & {F \in R^{n \times {({n - 1})}}} \end{bmatrix}}};$ the following is obtained from Q and Q⁻¹: ${{y({kp})}\overset{\bigtriangleup}{=}{\begin{bmatrix} {y_{1}({kp})} \\ {y_{2}({kp})} \\  \vdots \\ {y_{n}({kp})} \end{bmatrix} = {{Qx}({kp})}}};$ thus, y(kp+p)=Qx(kp+p), y(kp−p)=Qx(kp−p) and y(kp)=Qx(kp) are obtained, and the following is further obtained: ${\begin{bmatrix} {y\left( {{kp} + p} \right)} \\ {y({kp})} \end{bmatrix} = {\begin{bmatrix} {I_{n} - {\left( {p - \tau} \right)H}} & {{- \tau}H} \\ I_{n} & 0_{n} \end{bmatrix}\begin{bmatrix} {y({kp})} \\ {y\left( {{kp} - p} \right)} \end{bmatrix}}},$ where ${H = {{Q\overset{\_}{L}Q^{- 1}} = \begin{bmatrix} 0 & {C\overset{\_}{L}F} \\ 0_{n - 1} & {E\overset{\_}{L}F} \end{bmatrix}}};$ thus, the system is divided into two subsystems: $\begin{bmatrix} {y_{1}\left( {{kp} + p} \right)} \\ {y_{1}({kp})} \end{bmatrix} = {{\begin{bmatrix} 1 & 0 \\ 1 & 0 \end{bmatrix}\begin{bmatrix} {y_{1}({kp})} \\ {y_{1}\left( {{kp} - p} \right)} \end{bmatrix}} - {\begin{bmatrix} {\left( {p - \tau} \right)C\overset{\_}{L}F} & {\tau C\overset{\_}{L}F} \\ 0 & 0 \end{bmatrix}\begin{bmatrix} {\hat{y}({kp})} \\ {\hat{y}\left( {{kp} - p} \right)} \end{bmatrix}}}$ and ${\begin{bmatrix} {y_{1}\left( {{kp} + p} \right)} \\ {y_{1}({kp})} \end{bmatrix} = {\begin{bmatrix} {I_{n - 1} - {\left( {p - \tau} \right)E\overset{\_}{L}F}} & {{- \tau}E\overset{\_}{L}F} \\ I_{n - 1} & 0_{n - 1} \end{bmatrix}\begin{bmatrix} {\hat{y}\left( {{kp} + p} \right)} \\ {\hat{y}\left( {{kp} - p} \right)} \end{bmatrix}}},$ where ŷ(kp)=[y ₂(kp),y ₃(kp), . . . ,y _(n)(kp)]^(T) and ŷ(kp−p)=[y ₂(kp−p),y ₃(kp−p), . . . ,y _(n)(kp−p)]^(T);
 4. The distributed control method for consensus of an agricultural multi-agent system based on sampled data according to claim 3, wherein in step 3, the constraint conditions on the time delay and the sampling period that the multi-agent system achieves stability, that is, the sufficient and necessary conditions that the agents in the multi-agent system achieve average state consensus, are obtained by means of bilinearity and Hurwitz stability criteria, specifically: the following is obtained by using an invertible matrix T: ${{T^{- 1}E\overset{\_}{L}{FT}} = {\Lambda = \begin{bmatrix} \lambda_{2} & * & & \\  & \lambda_{3} & \ddots & \\  & & \ddots & * \\  & & & \lambda_{n} \end{bmatrix}}},$ where λ₂, λ₃, . . . , λ_(n) are non-zero eigenvalues of L and elements, * is 0 or 1, and further, {tilde over (y)}(kp+p)=T ⁻¹ ŷ(kp+p) {tilde over (y)}(kp)=T ⁻¹ ŷ(kp) {tilde over (y)}(kp−p)=T ⁻¹ ŷ(kp−p); a dimension reduction system is transformed into: ${\begin{bmatrix} {\overset{\sim}{y}\left( {{kp} + p} \right)} \\ {\overset{\sim}{y}({kp})} \end{bmatrix} = {\zeta\begin{bmatrix} {\overset{\sim}{y}\left( {{kp} + p} \right)} \\ {\overset{\sim}{y}\left( {{kp} - p} \right)} \end{bmatrix}}},$ where $\zeta = {\begin{bmatrix} {I_{n - 1} - {\left( {p - \tau} \right)\Lambda}} & {- {\tau\Lambda}} \\ I_{n - 1} & 0_{n - 1} \end{bmatrix}.}$ a characteristic polynomial of ζ is further obtained: ${{\det\left( {{sI}_{{2n} - 2} - \zeta} \right)} = {{❘\begin{matrix} {{sI}_{n - 1} - \left\lbrack {I_{n - 1} - {\left( {p - \tau} \right)\Lambda}} \right\rbrack} & {\tau\Lambda} \\ {- I_{n - 1}} & {sI}_{n - 1} \end{matrix}❘} = {❘\xi ❘}}};$ the following is obtained from the properties of a block matrix: ${{❘\xi ❘} = {{{❘{sI}_{n - 1}❘}{❘{{sI}_{n - 1} - \left\lbrack {I_{n - 1} - {\left( {p - \tau} \right)\Lambda}} \right\rbrack + {{I_{n - 1}\left( {sI}_{n - 1} \right)}^{- 1}{\tau\Lambda}}}❘}} = {{{\prod}_{i = 2}^{n}\left\lbrack {s^{2} - {\left( {1 - {p\lambda}_{l} + {\tau\lambda}_{l}} \right)s} + {\tau\lambda}_{l}} \right\rbrack} = {{\prod}_{i = 2}^{n}{g_{i}(s)}}}}};$ then, the following can be known by bilinear transformation of $s = {\frac{z + 1}{z - 1}:}$ f _(i)(z)=pλ _(i) z ²+2(1−τλ_(i))z+(p−2τ)λ_(i)+2.
 5. The distributed control method for consensus of an agricultural multi-agent system based on sampled data according to claim 4, wherein if f_(i)(z) is Hurwitz-stable and g_(i)(s) is Schur-stable, the stability is determined as follows: assuming z=ωι, then: f _(i)(ω)=−pλ _(i)ω²+2(1−τλ_(i))ωι+(p−2τ)λ_(i)+2, its real part is: f _(ω)(ω)=−pλ _(i)ω²+(p−2τ)λ_(i)+2 and its imaginary part is: f _(i)(ω)=2(1−τλ_(i))ω; upon verification, f_(i)(z) is Hurwitz-stable when satisfying the following conditions: ${\tau < \frac{1}{\lambda_{\max}}},{0 < p < {\frac{2}{\lambda_{\max}} + {2{\tau.}}}}$ that is, under the Hurwitz-stable condition, the multi-agent system achieves consensus.
 6. The distributed control method for consensus of an agricultural multi-agent system based on sampled data according to claim 5, wherein step 4 comprises the following steps: when L has an eigenvalue of 0, determining ϕ to have corresponding determinable eigenvalues of 0 and 1; the condition of Hurwitz stability of the system is determined after the tree transformation; when the constraint conditions are satisfied, the modulus values of remaining unknown eigenvalues of ϕ are less than 1, that is, the eigenvalues of ϕ except 0 and 1 are within a unit circle under the conditions; because G is undirected strongly connected, so a right eigenvector ω_(r) and a left eigenvector ω_(l) of ϕ when the eigenvalue is 1 satisfy the following: ${w_{r} = {w_{l} = {\frac{1}{\sqrt{n}}1_{n}}}},$ and ω_(l) ^(T) ω_(r)=1 is satisfied; all other eigenvalues of ϕ are within the unit circle and there is: ${\begin{bmatrix} {x\left( {{kp} + p} \right)} \\ {x({kp})} \end{bmatrix} = {{\phi\begin{bmatrix} {x({kp})} \\ {x\left( {{kp} - p} \right)} \end{bmatrix}} = {\phi^{k}\begin{bmatrix} {x(0)} \\ {x(0)} \end{bmatrix}}}},$ so ${{\lim\limits_{k\rightarrow\infty}\begin{bmatrix} {x\left( {{kp} + p} \right)} \\ {x({kp})} \end{bmatrix}} = {{\begin{bmatrix} {w_{r}w_{l}^{T}} & 0 \\ {w_{r}w_{l}^{T}} & 0 \end{bmatrix}\begin{bmatrix} {x(0)} \\ {x(0)} \end{bmatrix}} = \begin{bmatrix} {\frac{1}{n}{\sum}_{i = 1}^{n}{x_{i}(0)}} \\ {\frac{1}{n}{\sum}_{i = 1}^{n}x_{i}(0)} \end{bmatrix}}},$ therefore, the agricultural multi-agent system with a first-order dynamic model achieves average consensus despite time delays. 